Interplay of magnetic and structural transitions in Fe-based pnictide superconductors 
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The interplay between the structural and magnetic phase transitions occurring in the Fe-based 
pnictide superconductors is studied within a Ginzburg-Landau approach. We show that the magne- 
toelastic coupling between the corresponding order parameters is behind the salient features observed 
in the phase diagram of these systems. This naturally explains the coincidence of transition tem- 
peratures observed in some cases as well as the character (first or second-order) of the transitions. 
We also show that magnetoelastic coupling is the key ingredient determining the collinearity of the 
magnetic ordering, and we propose an experimental criterion to distinguish between a pure elastic 
from a spin-nematic-driven structural transition. 

PACS numbers: 74.70.Xa, 74.90. +n, 75.80.+q 



Introduction- The discovery of an unconventional 
high-temperature superconductivity in the F-doped 
arsenic-oxide LaFeAsO 1 has given rise to a great interest 
on iron pnictide compounds. To date several families of 
Fe-based superconductors have been discovered, and the 
superconducting transition temperature has been raised 
to above 50 K. 2 All of these systems display an intrigu- 
ing competition between structural, magnetic and super- 
conducting transitions. 2,3 The parent compounds of the 
so-called 1111 and 111 families undergo a structural tran- 
sition (ST) followed by a magnetic transition (MT) at a 
lower temperature, 4 whereas in the 122 and 11 cases these 
two transitions take place simultaneously. 5-7 In any case 
these orderings are quickly suppressed by doping or by 
applying pressure, which eventually gives rise to super- 
conductivity. The role, if any, of magnetic and elastic 
degrees of freedom in inducing this superconductivity is 
currently an open question. There is already a growing 
body of theoretical works advocating for spin fluctua- 
tion mediated superconductivity, 8 but the isotope effect 
observed in both magnetism and superconductivity sug- 
gests that the elastic medium also plays a role. 9 It is 
therefore compelling to understand the connection be- 
tween the ST and the MT in these systems. 

The ST, for example, has been studied in Refs. 10,11 
where it has been identified with a spin-nematic ordering, 
but ignoring the possible softening of the lattice itself. 
This idea has been further elaborated in Rcf. 12, where 
the observed softening of the lattice is interpreted as due 
to the fluctuations of the emerging nematic degrees of 
freedom. Recent studies 13 have also associated the ST 
with a possible ordering of the d xy and d yz Fe-orbitals in 
the Fe-As planes, with the consequent distortions having 
an impact on both the magnetic order 13 ' 14 and the elec- 
tronic structure. 15 On the other hand, the study of the 
MT in Ref. 16 focuses on the role of the magnetoelas- 
tic (ME) couplings in producing a weak first-order MT 
by means of the Larkin-Pikin mechanism. 17 This latter 
is relevant only in the fluctuation dominated Ginzburg 



regime and is subjected to the condition of a MT with 
diverging specific heat in the absence of ME coupling. 
Furthermore, the collinearity of the magnetic moments is 
examined in Ref. 18 by using a purely electronic model. 

In this paper we study the ST and the MT using a 
Ginzburg-Landau approach, in which the interplay be- 
tween elastic and magnetic degrees of freedom is consid- 
ered explicitly. This provides a general unified framework 
that goes beyond previous phenomenological approaches 
and rationalizes different experimental findings. We ad- 
dress, in particular, the simultaneity of the ST and MT, 
their character (first versus second order) , the collinearity 
of the magnetic structure, and the spin-nematic scenario 
for the ST. Our study identifies a particular ME coupling 
[see Eq. (2) below] as the common key factor behind the 
salient features of the ST and MT in the Fe-based super- 
conductors. 

Our main results are the following, (i) We derive the 
general phase diagram of the ST and MT, which exhibits 
four qualitatively different regimes, as shown in Fig. 1. 
This explains why the simultaneous ST and MT in the 
122 compounds are sometimes observed as second-order 
transitions 6 (la in Fig. 1) and sometimes as first-order 5 
(lb in Fig. 1). We moreover predict a first-order MT 
when this occurs separate but in the immediate vicinity 
of the ST (Ha in Fig. 1). This richness of the phase 
diagram is due to the ME coupling, (ii) The collinearity 
of the magnetic moments is linked to the existence of a 
particular ME term allowed by symmetry, (iii) An ex- 
perimental criterion is derived to distinguish between a 
pure elastic from a spin-nematic-driven ST. (iv) The fluc- 
tuations associated with ST are shown to become critical 
only along certain lines of high symmetry in the Brillouin 
zone. Along these lines the sound velocity vanishes in a 
pure elastic ST or the spin-nematic excitations are mass- 
less in a spin-nematic-driven ST. In both scenarios the 
ST has a mean-field behavior with no divergences. 11 

Ginzburg-Landau free energy - The basic ingredients to 
build up a Ginzburg-Landau theory for the ST and the 
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FIG. 1: Schematic T vs. doping diagrams of the possible 
structural and magnetic phase transitions in Fe-based su- 
perconductors (top panels) and general phase diagram T vs. 
Tg - T% determined by our Ginzburg-Landau theory Eq. (1) 
(bottom panel) . Empirically Tg — T^r increases with doping. 
In the regions I the structural and magnetic are simultaneous, 
while in regions II they take place at different transition tem- 
peratures, Tg and Tjv respectively. Continuous (dotted) lines 
indicate second (first) order phase transitions. The strength 
of the first order discontinuity is indicated by the intensity of 
the background color. 

MT in the Fe-based superconductors are the following. 
The initial order parameter for the ST can be taken as 
the shear strain u xy , 19 where subscripts refer to the prin- 
cipal axes of the tetragonal lattice (with two and four 
Fe atoms per unit cell for the 1111 and the 122 sys- 
tems respectively, see Fig. 2a). Thus u xy ^ implies 
the monoclinic distortion of the tetragonal unit cell ob- 
served experimentally, as illustrated in Fig. 2b. For the 
case of a spin-nematic driven ST an additional order pa- 
rameter a is required, which will be explained in more 
detail below. On the other hand the description of the 
MT requires two Neel vectors, 20 Li and L2, which can 
be associated with the magnetizations Mi and M2 of 
the two inter-penetrating Fe sublattices (Fig. 2b). The 
Ginzburg-Landau free energy then can be conveniently 
written as 

Fgl = Fm + Fe + Fme- (1) 

Here the magnetic part is F M = \A(L\ + L%) + \B(L\ + 
L|). In principle fourth order terms of the form L\L\ 
and (Li • L 2 ) 2 are also allowed by symmetry, but we 
omit them because the former does not affect the re- 
sults qualitatively while the latter is generated by the 
ME coupling (see below). For the elastic part we con- 
sider explicitly the critical strain u xy only and we write 
Fe = \c§QU xy + \f3u xy - The coefficients A and cqq are as- 
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FIG. 2: Tetragonal Fe-lattice (drawn as two interpenetrating 
sublattices) at temperatures T > Ts,Tn (a). A monoclinic 
distortion of the tetragonal unit cell and a collinear Neel order 
in the two Fe sublattices take place at T < Ts,Tn (b). 

sumed to vary with the temperature as A = A'(T — T^) 
and C66 = Cg 6 (T — Tg). and Tp. are the nominal 
MT and ST temperatures respectively (not to be con- 
fused with the actual transition temperatures Tjy and 
Ts), taken as the control parameters of our theory. It is 
our empirical obsevation that the results of the Landau- 
Ginzburg theory are consistent with the current experi- 
ments if we assume that T<1 — increses with doping. 
A possible expalantion of this trend at very low doping is 
the following. It is likely that the effect of doping is more 
on the magnetic sector (say, via the loss of nesting be- 
tween electron and hole bands) than the elastic medium 
which it affects indirectly. Thus, it may not be unrea- 
sonable to expect that at low doping Tg stays relatively 
unchanged while decreases (since the system hardens 
magnetically). The remaining coefficients B, C and /3 are 
taken as positive constants. As regards the ME term, the 
key contribution is 11 ' 12 ' 16 

Fme = 9iu xy {l^i ■ L 2 ), (2) 

where g% > in order to be consistent with the ex- 
perimental observation that the ferromagnetic Fe-Fe 
bonds are shorter than the antiferromagnetic ones in 
the collinear-Neel state. To our purposes, the standard 
magneto-striction uu(L\ + L|) can be neglected because 
it does not change qualitatively the results. It is worth 
mentioning that, in the context of the cuprate supercon- 
ductors, a similar free energy Fgl was used to study the 
structural transition of doped L^SrCuCv 21 

Phase diagram - We are in a position to discuss the 
general phase diagram for the ST and MT in the Fe- 
based superconductors (Fig. 1). We anticipate that, ex- 
perimentally, undoped and lightly doped 122 systems fall 
in regimes la and lb with simultaneous transitions of sec- 
ond and first order respectively. On the other hand, the 
1111 and 111 compounds appear to fall in the regime lib, 
where the transitions are separate (Tg > Tjy) and second 
order. In the regime Ila the ST is second order while the 
separate MT is first order. This regime has not yet been 
reported experimentally, perhaps because it is narrow in 
the tuning of physical parameters (such as pressure and 
doping) and therefore difficult to detect. 

Within our framework, the fact that these transitions 
sometimes are observed separate and sometimes simulta- 
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neous is explained as follows. Note that, once the mag- 
netic order sets in and the collinear-Neel state is formed, 
there is an effective shear stress that produces a mon- 
oclinic distortion via the ME term (Eq. 2). Thus, the 
MT in these systems implies a ST, while the converse 
is evidently not true. This gives rise to the regimes I 
and II in the phase diagram shown in Fig. 1, which are 
defined by the conditions < T$ and T§ > 1% respec- 
tively. Thus, as the temperature is lowered in regime II 
there is first the ST at T — Tg, where cqq = 0, followed 
by a separate MT at T N = T% + [giu (T N )]/A' , where 
Uo(T) = ±[(T| - TJcgg//?] 1 / 2 is the value of the strain 
order parameter u xy in the monoclinic phase at T < Tg. 
In the regime I, however, the system encounters the MT 
first at T = Tjy, where A = 0, but simultaneously the ST 
occurs due to the effective stress in the ME term. In this 
regime, unlike in II, there is no instability of the lattice 
because the elastic modulus cqq stays finite (and jumps 
across the transition). 

These two regimes are divided in their turn into sub- 
regimes la, lb, Ha and lib according to the character, 
first or second order, of the MT. This is again due to 
the ME coupling, and has nothing to do with symme- 
try reasons (such as the existence of cubic invariants). 
The point is that the effective coefficient of fourth-order 
term in the magnetic sector can become negative due 
to the ME coupling (in which case, one has to include 
higher order terms to bound the free energy from be- 
low). This is readily seen in the regime I (Tg < T$-), 
where the elastic sector is stable and one can ignore 
the u xy term in Fe- Indeed minimizing Fql one gets 
u xy = — <7i(Li • Li2)/c66, which implies that the elas- 
tic sector generates an effective magnetic term of the 
form -gi(Li ■ L 2 ) 2 /(2c 66 ). The MT into the collinear- 
Neel state turns first order when this term is sufficiently 
strong, with the tri-critical point given by cqq = g\/B 
and A = 0. The logic is quite similar in the regime II 
(Tg > Tjy), where one first has to expand u xy = uq + Suq 
in the monoclinic phase to then obtain (after further 
minimization) Suo = — 5i(Lj • L 2 )/(2|cg@|). This implies 
an effective magnetic term — gf(Li ■ L 2 ) 2 /(4|c 66 |) which 
drives the MT first order, with the tri-critical point given 
by A = giUo(Tv) and cqq = —g\/{2B). These results 
have to be contrasted to Ref. 16, where strong fluc- 
tuations are required and weaker first-order transitions, 
with narrower hysteresis, are expected. At the mean field 
level, it is the vicinity of the structural instability (where 
c 66 = 0) to determine the area (region Ha in the phase 
diagram of Fig. 1) where the MT becomes first order 
due to the ME coupling. In fact the discontinuity is the 
stronger the closer is to the point Tg = (as indicated 
by the color intensity in Fig. 1). 

Collinearity of the magnetic order. — At this stage, 
the role of the ME coupling in favoring the collinear- 
Neel state observed experimentally 22 is quite straightfor- 
ward. In regime II the collinearity is due to the term 
giu (Li ■ L 2 ), which lifts the degeneracy of the transi- 
tion temperatures for L± = -4= (Li ± L 2 ). This implies 



that only one between L± is non-zero in the ordered state 
(which one in particular depends on the sign of uo, i.e., on 
the ferroelastic domain). In regime I, however, the reason 
is different. In this case the set up of the collinear order 
is associated with the negative coefficient — g\j(2c§§) in 
front of the term (Li • L 2 ) 2 , generated by the ME cou- 
pling. Presumably, the vicinity to the structural instabil- 
ity makes this ME term strong enough to overcome any 
other contribution (due to, e.g., magnetic interactions) 
and thus forces the collinearity. 

Spin nematic versus pure elastic instability. — In the 
region II of the phase diagram, where the two transitions 
are separate, there are two different physical scenarios 
to conceptualize the ST: a pure ferroelastic transition or 
a transition induced by a spin-nematic ordering. In the 
following we discuss some physical implications that may 
allow to determine experimentally which one of the two 
scenario is realized in the Fe-based superconductors. 

(a) In the first scenario we have a proper ferroelastic in- 
stability in Fe caused by the vanishing of the elastic mod- 
ulus C66 (as described above) . An immediate consequence 
is that the temperature range over which the monoclinic 
distortion has square-root behavior uq(T) oc (Tg — T) 1 / 2 
is the same as the temperature range over which the 
elastic modulus is expected to have a linear dependence 
c 66 oc \T-T S \. 

(b) The second scenario is more subtle, as it arises from 
the possibility that, above Tjv, the system enters a spin 
nematic state in which the two Neel vectors are zero in 
average but fluctuate in phase, i.e., a = (Li • L 2 ) ^ 0. 23 
In this case the system is expected to undergo a pseudo- 
proper ferroelastic transition driven by the nematic or- 
dering. This scenario can be described by the free energy 

? a 2 , & 4 , C 66 2 i P 4 , (o\ 

F G L = + + ~Y U xy + 1 U xy + 92U xy (X, (3) 

where a = a'(T — T®) with T° being the nominal nematic 
transition temperature, while all other coefficients are 
T-indcpendent (including cqq for this discussion). The 
structural-nematic transition, below which u xy oc a 7^ 0, 
occurs at Tg = T® + 5|/( c 66 a where the effective elastic 
modulus cg6 = CQQ—g^/a vanishes. In this case, the range 
of temperatures over which cqq oc |T — Tg\ is restricted 
to the condition \T — Tg\ < g^Kc^a'). This range is in 
principle different from the T-range over which the mon- 
oclinic distortion has square-root dependence, as we have 
discussed in point (a) above. 

This simple observation, which is schematically rep- 
resent in Fig. 3, can be used as a quantitative cri- 
terion to distinguish between the pure ferroelastic and 
spin-nematic-induced scenarios. These considerations 
are based on a mean field analysis of the ST, which we 
justify in the following. 

Mean-field behavior for the structural transition - Pre- 
vious studies on ferroelastics 24 have shown that a second- 
order ST, like in the proper ferroelastic scenario (a), 
displays a mean-field behavior. This follows from the 
fact that at the transition point, where the elastic mod- 
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FIG. 3: Expected dependence on temperature of the squared 
spontaneous strain u 2 H and the elastic modulus C66 in scenar- 
ios (a) proper ferroelastic instability and (b) driven by spin- 
nematic ordering. Large deviations from a linear behavior in 
C66 as compared to u^. y indicate spin-nematic ordering. 

ulus c 66 = 0, the phonon velocity remains finite ev- 
erywhere in the Brillouin zone (q-space), except along 
the two lines of high symmetry q x = and q v = 
on the q z = plane. Consequently, except for these 
"soft" lines, the long-wavelength critical excitations as- 
sociated with the strain field u xy are gapped and the 
transition in three dimensions is mean- field like. 24 The 
contribution of these excitations to Fe can be written as 
SFe = \ E q ^o c 66(q)K a (q)| 2 + • • • , where the ellipsis 
denotes non-critical modes and interaction terms, 

C66(q) ~ C66 + ci cos 2 9 + C2 sin 2 2ip sin 4 9 + Kq 2 (4) 

near the soft lines, and 9 and (p represent the polar and 
the azimuthal angles respectively in q-space. The nq 2 
contribution is due to harmonic terms in the free energy 
with higher order derivatives (the usual stiffness associ- 
ated with collective excitations) and c\ and C2 are con- 
stants that depend on the elastic moduli of the system. 25 
In the scenario (b) the situation is slightly more sub- 
tle. Indeed, if we ignore the coupling g2U xy a, it appears 
that the nematic transition, which belongs to the Ising 
class, is in three dimensions below its upper critical di- 
mension and therefore it is not mean- field type. How- 
ever, the ME coupling gives rise to anisotropic correc- 



tions to the mass of the er-field. In fact the resulting 
mass term can be written as F a — \ ^^ q _^ a (q) | cr (q) | 2 , 
where a(q) — a — g 2 /cee(q). As a consequence, in this 
scenario too, the long-wavelength critical excitations are 
gapped except along two soft lines. 26 So the ST is com- 
pletely mean field in any case. 

We finally discuss some limitations of our study. We 
have restricted ourselves to magnetic states describable 
in terms of the Neel vectors Li and L2. This includes the 
collincar (it, 0) [or (0, it)] "stripe" order observed experi- 
mentally. In the case of a different magnetic ground state 
the corresponding Ginzburg-Landau theory may have a 
different form, with structural and magnetic orderings 
not necessarily coupled. This seems to be the case in 
BaMn2As2, where a G-type antiferromagnetic order is 
found and no structural distortion is observed. 27 The 
same applies to 11 Fe-chalcogenides, for which, in ad- 
dition, the question of whether the magnetic order is in- 
trinsic or arises from non-stoichiometry, i.e., due to inter- 
stitial Fe-ions inducing weak charge localization, is still 
open. 28 Critical fluctuations are also not considered in 
our theory. While we expect that they will not change 
the nature of the ST, they certainly will enhance the ef- 
fects of the Larkin-Pikin mechanism 17 on the MT. We 
do expect therefore an enlargement of the first order MT 
region (lb and Ha in Fig. 1). All these considerations 
deserve deeper investigation in future studies. 

In summary, we have shown that the magneto-elastic 
coupling Eq. (2) plays a key role in producing the spe- 
cific features of the structural and magnetic transitions 
in the Fe-based superconductors. This coupling is re- 
sponsible for the simultaneity of the transitions observed 
in the 122 systems and explains the different characters 
(first-order and second-order) of the transitions reported 
so far. It also naturally explains the collinearity of the 
magnetic order. We have also addressed the question of 
the possibility of a spin-nematic driven structural transi- 
tion, indicating an experimental way to discriminate this 
type of scenario from a proper ferroelastic transition. 

We acknowledge G. Garbarino, E. Kats, M. Nunez- 
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cussions. 
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